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Abstract 

We review our ongoing theoretical and phenomenological study of the discovery potential 
for instanton-induced DIS events at HERA. Constraints from recent lattice simulations will 
be exploited and translated into a "fiducial" kinematical region for our predictions of the 
instanton-induced DIS cross-section. 



*Plenary talk at the 3rd UK Phenomenology Workshop on HERA Physics, Sep 20-25, 1998, Durham/UK. 



1 Instantons 



Non-abelian gauge theories like QCD are known to exhibit a rich vacuum structure. The latter 
includes topologically non-trivial fluctuations of the gauge fields, carrying an integer topological 
charge Q. The simplest building blocks of topological structure are instantons (Q = +1) and 
anti- instantons (Q = —1) which are well-known explicit solutions of the euclidean field equations 
in four dimensions 0. 

Instantons (J) are widely believed to play an important role in various long-distance aspects @ 
of QCD: 

First of all, they may provide a solution of the famous Ua(X) problem (m r/ / 3> m v ), with the 
corresponding pseudoscalar mass splitting related to the topological susceptibility in the pure 
gauge theory by the well-known Witten-Veneziano formula ||. Moreover, a number of authors 
have attributed a strong connection of instantons with chiral symmetry breaking |||2|] as well as 
the hadron and glueball spectrum. 

However, there are also very important short-distance implications |]6|j7|j8|,p|,p!0[| of QCD instantons 
to which the present report is devoted: 

Instantons are known to induce certain processes which violate chirality in accord with the gen- 
eral axial-anomaly relation || and which are forbidden in conventional perturbation theory. Of 
particular interest in this context is the deep inelastic scattering (DIS) regime. Here, hard instan- 
ton-induced processes may both be calculated P,P,|10| within instanton-perturbation theory and 
possibly be detected experimentally []7| JTT] , [T^| , P~5[1 . As a key feature it has recently been shown [Q, 
that in deep-inelastic scattering (DIS) the generic hard scale Q cuts off instantons with large size 
p ^> over which one has no control theoretically. 

Our finalized results for inclusive instanton-induced DIS cross-sections are summarized in 



sections 2 and 4. Their weak residual renormalization-scale dependence is quite remarkable. 

As a second main point of this review (section 3), constraints from recent lattice simulations will 
be exploited |9,i~4jl and translated into a "fiducial" kinematical region for our predictions of the 



instanton-induced DIS cross-section based on instanton-perturbation theory. In section 5 we dis- 



cuss the expected event signature and search strategies based on our Monte Carlo generator 
QCDINS 1.60. Finally (section 6), we briefly address an interesting class of "fireball" events, ob- 
served in photoproduction, in the context of instantons and put forward a promising proposal |14| 



on extending our theoretical predictions beyond the regime of strict instanton perturbation theory. 



2 DIS cross-sections in instanton-perturbation theory 

In /-perturbation theory one expands the relevant Green's functions about the known, classical 
instanton solution = Ajp + . . . instead of the usual (trivial) field configuration = and 
obtains a corresponding set of modified Feynman rules. Like in conventional pQCD, the gauge 
coupling a s has to be small. 

The leading instanton-induced process in the DIS regime of e ± P scattering for large photon vir- 
tually Q 2 is illustrated in figure [l]. The dashed box emphasizes the so-called instanton- subprocess 
with its own Bjorken variables, 

Q' 2 = -q' 2 >0; x' = -^<l. (1) 
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Figure 1: The leading instanton- induced process in the DIS regime of e ± P scattering (rif = 3). 



It induces a total chirality violation A chirality = 2n/, in accord with the corresponding axial 
anomaly ||. In the Bjorken limit of /-perturbation theory, the dominant /-induced contribution 
to the inclusive HERA cross-section may be shown to take the form PJT0|] 
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The differential luminosity, d£^, accounting for the number of qg collisions per eP collision, has 
a convolution-like structure. It involves integrations over the gluon density, the 7-fl.ux P 7 and the 
known g-flux P^ 7 -* in the /-background (c.f. figure |]). The crucial instanton-dynamics resides in 
the /-subprocess total cross-section a^J(Q',x'), on which we focus our attention next pilO|]. 

Being an observable, a^(Q',x') involves integrations over all /(/)- "collective coordinates", in- 
cluding the / (Z)-sizes p (p) and the //-distance^ 4- vector R^, 
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The p(p)-integrals in (|3|) involve as generic weight the /(/)-density D(p(p)) [|3| JT5| . p~6 
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(4) 



with renormalization scale p r and N c = 3. 

^oth an instanton and an anti-instanton enter here, since cross sections result from taking the modulus squared 
of an amplitude in the single I background. 
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The function Q(R 2 / (pp) , . . .) in equation @, appearing in the exponent with a large numeri- 
cal coefficient 4ir/a s , incorporates the effects of final-state gluons. Within strict /-perturbation 
theory, it is given in form of a perturbative expansion fPT |, while in the so-called //-valley ap- 
proximation [18, 19 1 Q is associated with an analytically known closed expression [19|20| for the 
interaction between / and /, O ~ a s /(4:7r)S[A I IJ I ] — 1. With both methods agreeing for larger 
values of R 2 /(pp), we have actually used the valley method in our quantitative evaluation. 

Due to the nonvanishing virtuality Q' 2 in DIS, the "form factor" exp [—Q'(p +~p)] in (|3|), being 
associated with the off-shell quark (zero mode) q, suppresses large-size instantons Hence, 
the integrals in (|3|) are finite. In fact, they are dominated by a unique saddle-point PJlOll, 



p =p 



1/Q'; R* 2 ~l/(p + q 
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(5) 



from which it becomes apparent that the virtuality Q' controls the effective /-size, while x' 
determines the effective //-distance (in units of the size p). 

In figure 0, the resulting /-subprocess cross-sections (|3]) is displayed || over a large range of 
p r /Q' for fixed x' = 0.5 and Q'/A = 30,50,70. Apparently, we have achieved great progress 
in stability and hence predictivity by using the improved expression (^) of the /-density D(p), 
which is renormalization-group (RG) invariant at the 2-loop level, i.e. D~ l dD/d\n(p r ) = 0(a 2 ) . 
The residual dependence on the renormalization scale p r is remarkably flat and turns out to be 
strongly reduced as compared to the 1-loop case! Throughout, we choose as the "best scale", 
p r = 0.15 Q', for which da^/dp r ~ (c. f. figure ^). This choice agrees well with the intuitive 
expectation EE p r ~ l/(p) ~ Q' / f3 = O(0. 1) Q' . 
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Figure 2: Illustration of the weak residual renormalization-scale (p r ) dependence of the resulting 
/-subprocess cross-section a^J(Q',x'). 
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"Fiducial" region from lattice simulations 



There has been much recent activity in the lattice community to "measure" topological fluctu- 
ations in lattice simulations |2T] of QCD. Being independent of perturbation theory, such simu- 
lations provide "snapshots" of the QCD vacuum including all possible non-perturbative features 
like instantons (figure |3|). Let us discuss next, how these lattice results may be exploited to 
provide crucial support for the theoretical basis of our calculations in DIS: 

To this end, we first perform a quantitative confrontation P,T1| of the predictions from I- 



perturbation theory with a recent high-quality lattice simulation |23] of QCD (without fermions, 
rif = 0). The striking agreement which we shall find over a range of /-collective coordinates is a 
very interesting result by itself. 

Next, we recall (c.f. @ and (|5])) that the collective coordinate integrals in our DIS cross- 
section a^(Q',x') are dominated by a unique, calculable saddle-point (p*,R*/p*), in one-to-one 
correspondence to the conjugate momentum variables (Q', x'). This fact then allows us to translate 
the extracted range of validity of /-perturbation theory and the dilute /-gas approximation, 
(p < Pmax,-R/p > (-R/p)min), directly into a "fiducial" kinematical region {Q 1 > Q' min ,x' > x' min ) 
in momentum space! 

In lattice simulations 4d- Euclidean space-time is made discrete; specifically, the recent "data" from 
the UKQCD collaboration [23], which we shall use here, involve a lattice spacing a = 0.055 — 0.1 
fm and a volume V = ^ S p ace -/time = [16 3 -48 — 32 3 -64] a 4 . In principle, such a lattice allows to study 
the properties of an ensemble of /'s and /'s with sizes a < p < V 1 ^. However, in order to make 
instanton effects visible, a certain "cooling" procedure has to be applied first. It is designed to 
filter out (dominating) fluctuations of short wavelength 0(a) (c. f. figure |3] (a)), while affecting the 
topological fluctuations of much longer wavelength p ^> a comparatively little. After "cooling" , 
/'s and /'s can clearly be seen (and studied) as bumps in the lagrange density and the topological 
charge density (figure |] (b), (c)). For a more detailed discussion of lattice-specific caveats, like 
possible lattice artefacts and the dependence of results on "cooling" etc., see Refs. |2T|j23l . 

Of course, one has to extrapolate the lattice observables to the continuum (a =^> 0), before a 
meaningful comparison with /-perturbation theory can be made. This is complicated by a strong 
dependence of the various distributions on the number n coo i s of cooling sweeps for fixed (3 = 



Lagrange Density 



Topological Charge Density 




Figure 3: Instanton content of a typical slice of a gluon configuration at fixed x,y as a function of 
z and t [23|. (a) Lagrange density before "cooling", with fluctuations of short wavelength 0(a) 



dominating. After "cooling" by 25 steps, 3 /'s and 2 /'s may be clearly identified in the lagrange 
density (b) and the topological charge density (c). 
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In ref. p3] , however, equivalent pairs ((3,n coo \ s ) were found, for which shape and normalization 
of the distributions essentially remain invariant. For instance, the continuum extrapolation of 



the data for the (J + J)-density -D /+ j at n, 
be performed quite reliably [f 



'■cools I 



(6.0,23), (6.2,46), 



by simply rescaling the arguments p 



(6.4,80), may thus 
p(0)/p(a) ■ p. Here, 

_ A 



p(0) denotes the continuum limit of the weakly varying average p values, p(a), of D I+ j(p,a 
linear extrapolation in (a/Vo) 2 was employed. For consistency and minimization of uncertainties, 
one should use only a single dimensionful quantity to relate lattice units and physical units. 
Throughout our analysis, all dimensions are therefore expressed by the so-called Sommer scale p3, 



The resulting "continuum 
j ready to perform a 
For reasons of space, 



25 ro, with 2ro ~ 1 fm, which we prefer over the string tension |23 

data" for D I+ j(p) are displayed in figure f|. They scale nicely. We are now ready to perform a 
quantitative comparison with the predictions of /-perturbation theory (14 
let us concentrate here on the (/ + /)-density D I+ j(p). The prediction (|]) of /-perturbation theory 
is a power law for small p, i. e. approximately D ~ p 6 for nj = 0. Due to its 2- loop RG-invariance 
the normalization of D I+ j(p) is practically independent of the renormalization scale p r over a 
wide range. It is strongly and exclusively dependent on r A^ n ^ =0 , for which we take the most 
recent, accurate result by the ALPHA-collaboration |]2"5| , 2r Aj^g n ^. =0 = (238 ± 19) MeVfm. In 
figure ^ (b) we display both this parameter-free prediction from (|4j) of /-perturbation theory and 
the continuum limit of the UKQCD data in a log-log plot, to clearly exhibit the expected power 
law in p. The agreement in shape and normalization for p<0.3 (2rn.) — 0.3 fm is striking, indeed, 
notably in view of the often criticized "cooling" procedure and the strong sensitivity to Aj^g n ^ =0 . 



By a similar analysis []I4| |, we were able to infer from the "equivalent" UKQCD lattice data a range 
of validity -R/p> 1 of the valley expression for the //-interaction Q(R 2 /(pp), . . .) in (^). Finally, 
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Figure 4: Continuum limit [|E4j] of "equivalent" UKQCD data [23] for the (/ + /)-density at 
(/?, "cools) = (6.0,23) [□], (6.2,46) [o], (6.4,80) [A]. The striking agreement with 2D(p) of /- 
perturbation theory from (^) is apparent [l^|. The 3-loop form of a s with Aj^g n =0 from AL- 
PHA was used, (a) For p r = 1.2/ p, the agreement extends up to the peak; (b) Log-log plot 
to exhibit the expected power law ~ p 6 and the agreement in magnitude for small p over a wide 
range of p r . The dashed error band results from varying A^jg nj!=0 and p r within its error and 
given range, respectively. 
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we have confirmed p3| , |14f the approximate validity of the dilute-gas picture for sufficiently small 
instantons^ with p<(0.3 — 0.5) fm. The latter results are based on the "packing fraction" |2S 



being < 1 and a test of the dilute-gas identity: (Q 2 ) = N tot . Here Q is the topological charge and 
N tot the total number of charges. These results strongly support the reliability of our calculations 
in DIS. 

By means of the discussed saddle-point correspondence (||), these lattice constraints may be 
converted into a "fiducial" region for our cross-section predictions in DIS M, 
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(6) 



4 HERA cross-section 

Figure ^] displays our finalized /-induced cross-section at HERA P,|TU||, as function of the cuts 
x' mhl and Q' min . For the minimal cuts @ extracted from the UKQCD lattice simulation, we obtain 
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Figure 5: /-induced cross-section at HERA as function of the cuts in (V, Q'). 

a surprisingly large cross-section, 

^HERA<y > 0.35, Q' > 8 GeV) ~ 126 pb; x Bj > 1(T 3 ; 0.9 > y Bj > 0.1. (7) 

Hence, with the total luminosity accumulated by experiments at HERA, £ = O(80) pb -1 , one 
already expects O(10 4 ) /-induced events on tape from this kinematical region. Note also that 
the cross-section quoted in Eq. (^) corresponds to a fraction of /-induced to normal DIS (nDIS) 
events of 



f 



(I) ^HERA 
" (nDIS) 
°HERA 



0(1)%; 



for x Bj > W~ s ; 0.9 > s/Bj > 0.1. 



2 Note that the full (/ + 7)-ensemble without the size restriction is known not to be a dilute gas 
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This is remarkably close to the published upper limits on the fraction of /-induced events [26 
which are also on the one percent level. 

There are still a number of significant uncertainties in our result for the cross-section. For fixed 
Q' and x' cuts, one of the dominant uncertainties arises from the experimental uncertainty in the 
QCD scale A. In the 2- loop expression for a s with rif = 3 (massless) flavours we used the value 
= 282 MeV, corresponding to the central value of the DIS average for n f = 4, A^| = 234 

MeV [27]. If we change A^ within the allowed range, ~ ±65 MeV, the cross-section (0) varies 
between 26 pb and 426 pb. Minor uncertainties are associated with the residual renormalization- 
scale dependence (c.f. figure |2|) and the choice of the factorization scale. Upon varying the latter 
by an order of magnitude, the changes are in the O(20) % range only. 

By far the dominant uncertainty in o"j/j; RA arises, however, from the uncertainty in placing the 
(x',Q ; ) cuts (c.f. figure |5|). Hence, the constraints (|6]) from lattice simulations are extremely 
valuable for making concrete and reliable predictions of the /-induced rate at HERA. 



5 Signatures and searches 

An indispensable tool for investigating the structure of the /-induced final state and for develop- 
ping optimized search strategies is our Monte-Carlo generator for /-induced DIS-events, QCDINS 
1.60. Besides the matrix element for the /-induced hard subprocess, it provides leading- log parton 
showers and hadronization via its interface to HERWIG 5.9. 

The characteristic features of the /-induced final state are illustrated in figure ^ (a) displaying 
the lego plot of a typical event from QCDINS 1.60 (c.f. also figure [1]): 




Figure 6: (a) Lego plot of a typical instanton-induced event from QCDINS 1.60. 

(b) An interesting real "fireball" event in photoproduction from ZEUS |[(| with very large total 

E T and multiplicity. 

Besides a single (not very hard) current- quark jet, one expects an accompanying densely populated 
"hadronic band". For £Bjmin — 10~ 3 , say, it is centered around Tj ~ 2 and has a width of Arj ~ ±1. 
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The band directly reflects the isotropic production of an /-induced "fireball" of O(10) partons 
in the /-rest system. Both the total transverse energy (Et) — 15 GeV and the charged particle 
multiplicity (n c ) ~ 13 in the band are far higher than in normal DIS events. Finally, each /- 
induced event has to contain strangeness (and possibly also charm) such that the number of K 0, s 
amounts to ~ 2.2/event. 

Despite the high expected rate (|7[) of /-induced events at HERA, no single observable is known 
(yet) with sufficient nDIS rejection. Hence, a dedicated multi-observable analysis seems to be 
required. Neural network filters are being tried and exhibit a very good analyzing power if 
applied to >C(5) observables [pSj] . Strategies to produce "instanton-enriched" data samples and 



to reconstruct (Q' 2 , x') are under study and look quite promising f28|. Clearly, in all cases, a good 
understanding of the perturbative QCD background in the tails of the considered distributions is 
required. 



6 Going beyond instanton-perturbation theory 

A class of striking "fireball" events in photoproduction, with large total Et and large multiplicity, 
has been reported |35| at this meeting (see e.g. figure |5] (b)). While the quantitative analysis 
is still in an early stage, these events seem to exhibit all characteristics of /-induced events 
(c. f. figure H (a)). However, it appears that - unlike ordinary QCD perturbation theory - 
the hard photoproduction limit, Q 2 =^ 0, (-Er)jet large, is not within the reach of strict /- 
perturbation theory. The reason is that in the Q 2 =^ limit, we encounter a contribution 
to the photoproduction cross-section, which tends to diverge, if integrated over the /-size p. 
This IR divergence at large p is independent of the Et of the (current quark) jet and directly 
associated with the "bad" large-p behaviour of the perturbative expression (f|) for the /-density, 
D(p) ~ p 6 ~ 2 / 3 "/. In contrast, the actual form of D(p) (c.f. figure |j) is strongly peaked around 
p ~ 0.5 fm, and appears to vanish exponentially fast for larger p. The above "fireball" events 
and more generally, the ongoing /-searches at HERA, provide plenty of motivation for trying to 
extend our calculational framework beyond strict /-perturbation theory. We are thus led to make 
the following promising proposal in this direction: 

One may try and replace the most strongly varying entries in the perturbative calculations, the 
/-density D(p) and the //-interaction il(R 2 /(pp), . . .), in @, by their actual form as extracted 
from the recent non-perturbative lattice results [|i~4j1 . 

The /-rates in photoproduction and the /-contributions to further interesting observables may 
then be calculated, and the (Q' 2 , x') cuts in /-searches be considerably relaxed! Due to the strong 
peaking of D(p), only the region around p ~ 0.5 fm enters and the dilute-gas approximation may 
well continue to hold up to the peak (c.f. section 3). 
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